MINIMAL CUBATURE RULES AND POLYNOMIAL 
INTERPOLATION IN TWO VARIABLES 



YUAN XU 



Abstract. Minimal cubature rules of degree An — 1 for the weight functions 

Wc.,p,±i y) = 1^ + - y\^^+\a - - 

on [—1, 1]^ are constructed explicitly and are shown to be closed related to the 
Gaussian cubature rules in a domain bounded by two lines and a parabola. 
Lagrange interpolation polynomials on the nodes of these cubature rules are 
constructed and their Lebesgue constants are determined. 



1. Introduction 

Minimal cubature rules have the smallest number of nodes among all cubature 
rules of the same precision. Let be a non-negative weight function on a domain 
^2 C M^. For a positive integer s, a cubature rule of precision s with respect to W 
is a finite sum that satisfies 

N 

(1.1) / f{x,y)W{x,y)dxdy = y2Xkf{xk,yk), V/ e H^, 
•^^ k=i 

where denotes the space of polynomials of degree at most s in two variables, 
and there exists at least one function /* in H^^i such that the equation (1.1 ) does 
not hold. 

It is known that the number of nodes of a cubature rule necessarily satisfies 

(1.2) N > dimn^,_i = ^ , s = 2n - 1 or 2n - 2 

(cf. [121 mi)- A cubature rule of degree s with N attaining the lower bound in 



( 1.2 ) is called Gaussian. Unless quadrature rules in one variable, Gaussian cubature 
rules rarely exist. At the moment, they are known to exist only in two cases. The 
first case, discovered in |13j. is for a family of weight functions that includes, in 
particular, Wa,i3,±i2 defined by 

(1.3) W^.,^±i{u,v) = (l-u + w)"(l + u + tO'^(M^-4t;)±5 

on the domain fl ~ {{u^v) : 1 + u + w > 0, 1 — it + v > 0, > 4u}, bounded by two 
lines and a parabola. On the other hand, Gaussian cubature rules of degree 2ri — 1 
do not exist when W is centrally symmetric, that is, when W and its domain are 
both symmetric with respect to the origin: (— a;, —y) £ O whenever {x,y) G and 
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W{—x, ~y) — W{x^y). For centrally symmetric weight fmictions and s = 2n — 1, 
a stronger lower bound 10 for the number of nodes is given by 

„o n nin +1) n \ 

(1-4) jV>dimn^i+ y = ^ ^ + [2J ■ 

A cubature rule that attains this lower bound is necessarily minimal. There are, 
however, only a couple of examples for which this lower bound is attained for all n, 
most notable being the product Chebyshev weight functions on the square. 

In the present paper we shall show that the minimal cubature rules of degree 
An — 1 exist for a family of weight functions that includes, in particular, 

(1.5) W„,^,±i (x,2/) := \x + - y\'f^+\l - x')^^! - yYK 

on [—1, 1]^ and, furthermore, there is a connection between these minimal cubature 
rules and Gaussian cubature rules associated with the weight function /3 ± i ■ 



The weight functions (1.5) include the product Chebyshev weight functions (when 
a = (3 — zt^), for which the minimal cubature rules are known to exist and have 
been established in several different methods [H [TTJ [T^ . Our result shows that 
they can be deduced from the Gaussian cubature rules for W_ 1 _ 1 _|_ 1 on il. Giving 
the fact that so few minimal cubature rules are known explicitly, this connection is 
rather surprising. 

Cubature rules are closely related to interpolation by polynomials. Based on the 
nodes of a Gaussian cubature rule of degree 2n — 1, there is a unique Lagrange in- 
terpolation polynomial of degree ri — 1 which converges to / in norm as n — ?> cxd 



f |18|). On the nodes of the minimal cubature rule that attains (1.4), there is a 
unique Lagrange interpolation polynomial in an appropriate subspace of polyno- 
mials jl9j . Furthermore, the interpolation polynomials based on the nodes of the 
minimal cubature rules for the produce Chebyshev weight function yV_ 1 _ 1 _ 1 , 

studied in [2^, has the Lebesgue constant of order (logn)^ ^3J, which is the min- 
imal order of projection operators on [—1,1]^ ^Tj]. We shall discuss the Lagrange 
interpolations based on both the nodes of Gaussian cubature rules with respect to 



W^^p ±i ^^-d the nodes of minimal cubature rules for ( 1.5 ) in this paper. 

The paper is organized as follows. The next section is the preliminary, in which 
we recall basics on cubature rules and, in particular, the connection between cu- 
bature rules and interpolation polynomials, as well as basics on the orthogonal 
polynomials that will be needed in the paper. The Gaussian cubature rules for 
weight functions including W^a,/3,±i2 and minimal cubature rules for W^ p j-i are 
discussed in Sections 3 and 4, respectively. The interpolation polynomials based on 
the nodes of these cubature rules are treated in Sections 5 and 6, respectively. 



2. Preliminary and Background 

Minimal cubature rules are closely connected to orthogonal polynomials and 
to polynomial interpolation. We recall the connections in this section and state 
necessary definitions and properties of the weight functions and their orthogonal 
polynomials that will be needed later in the paper. 

2.1. Cubature, orthogonal polynomials and interpolation. Let be a non- 
negative weight function defined on a domain fl in that has all finite moments; 
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that is, J^x-{x2W{xi,X2)dxidx2 < oo for all j,k e Nq. Then orthogonal polyno- 
mials of two variables with respect to W exist. Let Vn{W) denote the space of 
orthogonal polynomials of degree n in two variables. Then 

dimVniW) = n+l. 

Assume that W is normalized so that W{xi,X2)dxidx2 = 1. A basis of Vn{W), 
denoted by {Pk,n : < /c < n}, is mutually orthogonal if 

Pk,n{xi,X2)Pj,nixi,X2)W{xi,X2)dxidx2 = hkSk^j, < k,j < n, 

where hk > and it is called orthonormal if /ife = 1 for < fc < n. The reproducing 
kernel Kn{W; ■, ■) of Ul in L^{W) is defined by 

Kr.iW;x,y)piy)W{y)dy^p{x), Vp G , 

in which x — {xi,X2) and y — (yi,j/2)- If Pk,n are orthonormal, then the reproduc- 
ing kernel Kn{W; •, •) of II^j in L^{W) is given by 

n m 

(2.1) Kr,{W; x,y)^Y.Yl Pk,MPk,m{y)- 

m=0 k=0 

Recall that a Gaussian cubature rule of degree 2n — l, as in (1.1), has dimll^.j 
nodes. These nodes are necessarily common zeros of orthogonal polynomials in 
VniW), that is, zeros of all polynomials in Vn{W) (cf. [3 [lOl E] ) . 

Theorem 2.1. Let n > 1. A Gaussian cubature rule of degree 2n — 1 exists if and 
only if its nodes are common zeros of orthogonal polynomials of degree n. Moreover, 
the weights Xk of the Gaussian cubature rule are given by 

Xk = {K^a-i{W; {xk,yk), {xk,yk))r^ , l<k<N. 

Unlike interpolation in one variable, polynomial interpolation in two variables 
may not exist for a set of distinct points. It does exists if the interpolation points 
are nodes of a Gaussian cubature rule [T^ . 

Theorem 2.2. Let N — dimn^_]^ and let {{xk,yk) ■ ^ < k < N} be the nodes 
of a Gaussian cubature rule of degree 2n — 1. Then there is a unique interpolation 
polynomial, Lnf , of degree n — 1 that satisfies 

f{Xk ) ~ f{xk,n,yk,n), I < k < N. 

Furthermore, this interpolation polynomial is given explicitly by 

N 

Lnf{x,y) ^^f{xk,yk)ik{x,y), tk{x,y) := XkKn-l{W■,(x,y),{xk,yk))■ 
k=l 

For centrally symmetric weight functions, we consider the minimal cubature 
rules whose number of nodes attains the lower bound in (1.4 1. The nodes of such 
a cubature rule are common zeros of a subspace of Vn{W) (tlQj, see also [l9]). 

Theorem 2.3. Let n > 1. The minimal cubature rule of degree 2n—l that attains 
the lower bound (1.4) exists if and only if its nodes are common zeros of [^^y^J + 1 
orthogonal polynomials of degree n. 
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Since the number of nodes of such mimmal cubature rules is iV = dimn^_]^ + [^J , 
the polynomial that interpolates at the nodes of the cubature rule needs to be from 
a polynomial subspace of n„ that has dimension N. An obvious candidate of this 
subspace is the hnear span of \ In, where I„ :— {Qk.n : fc = 0, 1, . . . , L^^J } 
denotes a set of orthonormal polynomials that vanish on the nodes of the minimal 
cubature rule. Let {Pk,n '■ ^ ^ k < [f J} be the orthonormal basis of the orthogonal 
complement of in Vn{W). Then Pj^ n G V„(iy) and none of P^, „ vanishes on all 
nodes of the cubature rule. We define a subspace 11* of 11^ by 

(2.2) n; n2_^ U span {Pfc,„ : 1 < fc < 



The weights \k,n of the minimal cubature rule in Theorem 2.3 are given in the 
lemma below. 



Lemma 2.4. Let Pk.n be as in (2.2 I. There exists a sequence of positive numbers 
{bk,n : 1 < < LfJ}' uniquely determined, such that the kernel K^{-, •) defined by 

LfJ 

(2.3) K:iW;x,y) = Kn-i{W;x,y) +J2bk,nPkA^)PkAy)^ 

fc=i 

where x — {xi,X2) and y = (2/1,2/2); satisfies 

(2.4) \k,n^[K{W;{xk,yk),ixk,yk))r\ l<k<N. 

This lemma was proved in [T5] and the coefficients were shown to be determined 
by the matrix [Cn{Pj,nPk,n)]] k=o EH, where {Po,n, ■ ■ ■ ,Pn,n} is an orthonormal 
basis of V„(Ty) and C„/ — "^^^i ^kfi^k^Vk) is the minimal cubature rule. The 
kernel K*{-, •) can also be used for the Lagrange interpolation polynomials based 
on the nodes of the minimal cubature rules, as stated in the following theorem |19j. 

Theorem 2.5. Let W be a central symmetric weight function. Let N — A\mIi^^_i-\- 
[^J and let {(xk,yk) ■ ^ l£ k < N} be the nodes of the minimal cubature rule of 
degree 2n — 1. Then there is a unique interpolation polynomial, Cnf , in H* that 
satisfies 

Cnf{xk,n,ykji) = f{xkM,yk,n), I <k < N. 

Furthermore, this interpolation polynomial is given explicitly by 

N 

^nf{x,y) ^^f{xk,yk)h{x,y), ik{x,y) -.^ \k,nKl{W]{x,y),{xk,yk)), 

k=l 

where Xu.n o-f^ the weights of the cubature rule given in (2.4 1. 

2.2. Weight functions and orthogonal polynomials. We define the weight 
functions for our Gaussian and minimal cubature rules. Our first weight function 
is defined on the domain 

-.^ {{u,v) : 1 + u + V > 0,1 - u + V > 0,u'^ > Av}, 

bounded by a parabola and two lines, as depicted in Figure 1. Let w be a nonneg- 
ative weight function defined on [—1, 1]. We define 

(2.5) Wy{u,v) -.^ bi^^^w{x)w{y){u'^ - 4:vy , {u,v)efl, 
where the variables {x,y) and {u,v) are related by 

(2.6) u — X + y, V — xy 
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and fe^,^ is the normalization constant such that W^{u,v)dudv = 1. In the case 
of the Jacobi weight function w = Wa,i3 defined by 

the weight function W-y is denoted by Wa,p.-f and it is given by 
(2.7) W„^p^y{u,v):^bo,,fi^y{\~u^v)'^{l + u + vf{u^-'iv)\ {u,v)£n, 
where a,(5,"f > —1, a + 7 + | > — 1 and /3 + 7 + ^ > — 1 and [T4( Lemma 6.1], 
o w ^ r(a + /3 + 7+|)r(a + /3 + 27 + 3) 

^ ■ ^ "'^'^ ■ 22"+2^*+47+2 r(a + i)r(/3 + i)r(7 + i)r(a + 7 + |)r(/3 + 7 + §) ' 

The weight function is related to w{x)w{y) by the relation 

(2.9) / f{u,v)W^{u,v)dudv = b^^^ f{x + y,xy)w{x)w{y)\x - y\'^^~^'^dxdy, 
Jn J A 

where :— {{x^y) : —1 < x < y < 1}. Since the integral in the right hand side has 
symmetric integrand in x, y, it is equal to half of the integral over [—1,1]^. In partic- 
ular, if Cw is the normalization constant of w on [—1, 1] so that J w{x)dx — 1, 
then the normalization constant of W_ i is given by 6^, _ i = 2c^ . 

The orthogonal polynomials with respect to Wa,p,-y were first studied by Koorn- 
winder in [6 and further studied in 7,8, 14J. They were applied to study cubature 
rules in [T3]. In the case of 7 = the orthogonal polynomials with respect to 
Wj. 1 can be given explicitly. Let p„ denote the orthogonal polynomial of degree n 
with respect to w. Then an orthonormal basis with respect to W_i is given by 

^2 10) P^"^^(u v) = +P«(y)Pfc(a;), < k < n, 

' \V2pnix)Pn{y), k^n, 

and an orthonormal basis with respect to Wi is given by 

(2.11) p(i)^^^,)^ Pn+li-)Pkiy)~Pn+l{y)M-) ^ o<k<n, 

X y 



both families are defined under the mapping (2.6 1. In the case of Wa.p^-y we denote 

the orthogonal polynomials by P^'^'^ ■ In particular, P^'^ '^^ are expressible by 
the Jacobi polynomials. 

Our second family of weight functions are defined on [—1,1]^ by 

(2.12) W-,{x,y):^W^{2xy,x^ + y^-l)\x^-y% (x, y) G [-1, 1]^, 
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where is the weight function in (2.5). In the case of Wa,i3,-y, it becomes 
(2.13) y) := 6a,/3,^4^|a; - y\'"+'\x + y\'P+\l - x'^il - y^)\ 



which includes ( 1.5 ). The is normalized if is because of the integral relation 

(2.14) I f{u,v)W^{u,v)dudv= I f{2xy,x^ + y^-l)W^{x,y)dxdy. 
Jn "'[-1,1]^ 

The orthogonal polynomials with respect to W-^ can be expressed in terms of or- 
thogonal polynomials with respect to ([H])- For this paper we will only need 
a basis for V2n{W~f), which consists of polynomials 

,,,,, iQk.Li^,y)--=Pt'!li'2^y^^' + y''^): o<fc<,^, 

(2.15 , ■ , / s , , 

2Qi'Li-, y) - y')P^^:l'll\2xy, x^ + y'-l), 0<k<n-l, 

where P^^^^'li are orthonormal polynomials with respect to the weight function 
(1 — 7i + w)(l + 7i + v)W^{u, v) and bi^'^^ is a normalization constant for the weight 
function. In the case of 'yVa,^,^, we denote the orthogonal polynomials by , 

in which case, Pj.'^^'^i = Pkn-i^^'^ (2.151. We will need explicit formulas for 
these polynomials when 7 = ^1/2, which we sum up in the following subsection. 
Further results on orthogonal polynomials with respect to can be found in [22] . 

2.3. Jacobi polynomials and orthogonal polynomials for >V^ p _i. The Ja- 

cobi polynomials are orthogonal with respect to Wa^p and they are given explicitly 
by a hypergeometric function as 

^) ^ ^^.^1 (~"' " ; : t ^ ^) = + • • ■ ' 

where zi"'^-* is the leading coefRcient. By 16, (4.21.6)], 

^ ^ " 2"n! "'^ 2«+/^+ir(a + l)r(/3 + l) 

The Jacobi polynomials satisfy the orthogonality conditions 

Ca,/3 j'^P!C'^\x)Pt^Hx)w^A^)dx = hl^^'^^Sn,m, 

where 

,(a,0) ^ {a + l)nW + l)n{a + f^ + n + l) 
^ ' " ■ n!(a + /3 + 2)„(a + /3 + 2n+l) ' 

The reproducing kernel of kik'^^ of the space of polynomials of degree at most n is 
given by, according to the Christoffel-Darboux formula, 

{2n + a + p + 2)[a + l)„(/3 + 1)„ 

Pj^f (x)Pi"-^^(;/) - P^^f (y)Pi"-^^(:.) 
The Gaussian quadrature of degree 2n — 1 for the Jacobi weight is given by 



(2.19) c^.p 



C f{x)w^^p{x)dx = ^t^^fi^k.n), V/ e n2„_i, 
•^-1 k=l 
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where Xi^n, ■ ■ ■ , Xn^n are the zeros of the Jacobi polynomial Pjf'^'' and 

We denote the orthonormal Jacobi polynomials by p^''^\ It follows readily that 
pl^''^\x) = {hi^'^^)~^ Pn^'^^x). The following lemma will be needed in Section 6. 



Lemma 2.6. For a,l3 > —1 and m > 0, define 

n 

(2.20) /J„, :=^a("'«(1 



4) 



(q + 1,;3 + 1) 
Pn-1 



(Xk) 



k=l 



Then hm — Ca,p / Ca+i^p+i for < m < n — 2, and 

4(l + a)(l + /3)(l + a + ^ + 2n) 



hn-l 



(2 + a + /3)(3 + a + I3){1 + a + /3 + n)' 

Proof. For < m < n — 2, we can apply Gaussian quadrature and use the orthonor- 
mality of Pm'^^'^^^'^ to conclude, since (1 — x'^)'Wa.i3{x) — Wa+i,/3+i{x), 

= c„,^ (1 - x^) \pi:+''^'+'\x)Y w^A^)dx = ''"-^ . 

J-l L -I Ca+l.fS+l 

For TO = ri — 1, we cannot apply the Gaussian quadrature of degree 2rt — 1 directly, 
since (1 — x^)[pm^^'^^^\x)]'^ has degree 2n. However, by [16, (4.5.5)], 

Xk)K.-i ^^'=^-2n + a + /3)(2n + a + /3 + l)^"-i ^'''^^ 



4n(n + 1) 



2n + a + P + l)(2n + a + /S + 2) 
and by the three-term relation satisfied by {Pi°''^''}n>o [HI (4.5.1)], 
{n+l){n+a+f3+l){2n+a+l3)Pi°l:f\xk) = -(n+a)(n+/3)(2n+a+/3+2)p/,"f ^(a;^). 
From these two equations it follows that 



(2.21) {1 - xl)PtJ,'^'^'\x,) 



4(n + a){n + (3) 



{2n + a + I3){n + a + 13 + 1) 



Denote the coefficient in front of P^^i\xk) in the above equation by Z?„. Then, 
by the Gaussian quadrature and the orthogonality of the Jacobi polynomials. 



k=l 



eV'^^'^i c.,w p}rj,'''^'>(x)p^r-'i'(-)^.Ax)dx 



Ca.Ji- 



,(a + l.f} + l) ^1 ^ 
'n-1 



'n-1 



■Pi-f''(2;) Wa.Ax)dx 



,{a,0) ,(a+l,/3+l) 

"n-l 'n-1 

,(q+1,;9+1) ia,l3 ' 
"n-l 'n-1 



which simplifies, by (2.16) and (2.17), to the stated result for h 



□ 
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In Section 6 we will need the explicit formula of orthogonal polynomials and 
reproducing kernels for the weight function yVa,i3,-i/2i which we rename as 

I _ |2a+l| I ^|2/3+l 

(2.22) Wo.A^,y) W„.^,_i/2(a:, y) = 2cl/ ^L=^—M^. 

In this case the orthonormal polynomials of even degree in (2.15) are given in terms 
of Jacobi polynomials, which are 

Proposition 2.7. Let a,l3 > —1. An orthonormal basis o/V2n(VVQ/3 _i) *s given 
by, for < k < n and < fc < n — 1, respectively, 

lQi"2n (cose, cos <j,) = pt^\cosi0 - <j,))p^t^\cOs{e + <P)) 

+ p^^)(cos(0 - (/.))P^«(COS(0 + 0)), 

2Qt£{cose,coscl^) = 7a,/3(^' - y^) [ptV'^+'^(cos(0 - 0))pi"+i'^+i)(cos(0 + 0)) 

+pi"+^'^+^)(cos(0 - 0))pt+^'^+^)(cos(0 + 

where iQ^^.tn ^''^-d- 2Q^n,2n "'^'5 multiplied by a/2/2 and^a.fi — Ca+i^p+i/{V2ca,p)- 

Denote the reproducing kernel of 11^ with respect to yVa,i3 by /C"''^(-, •). By 
Theorem 4.8], the kernel /C2n-i given explicitly by 

(2.23) Ktl,{x,y)^K-\{s,t) + d^^f {xl - xlM - yl)K+l'^+\s,t) 

+ S^f{x,+x2){yi+y2)K^f^'{s,t) 
+ di\fix^-x2)iyi-y2)K:+l'^is,t), 

where s = {2x^X2, xj+x^- 1), t = (2yij/2, + - 1), d^^f = (t+i,b+j/ci,i} and, 

with {xi,X2) — (cos0,cos02) and (yi,y2) — (cos cos ^2), 

(2.24) 

K^'Hs, t) \ [fc^'''(cos(0i - 02), C0S((/.1 - (/)2))fc^^^(cos(0i + 02), COs(0i + ^2)) 
+ fc,^"3(cos(0i - 02), COS(01 + 02))fc:'^(cos(0i + 02), COs(0i - ^2))] ■ 

For orthonormal basis of odd degrees and the reproducing kernels of even degrees, 
as well as other results on them, see |22j . 

3. Gaussian cubature rules 

In this section we consider Gaussian cubature rules for W+i on and their 
transformations. We shall show that these rules can be transformed into minimal 
cubature rules for y\>a,ii-\l2 on [—1,1]^ in the next section. The first proof that 
Gaussian cubature rules exist for Wj^ i was given in |13j via the structure matrices 
of orthogonal polynomials. Below is another proof that is of independent interest. 

We start with the Gaussian quadrature rule for the integral against w on [—1,1], 

/I " 
f{x)w{x)dx = ^ \kfixk,n), f e n2„_i, 

'1 k=l 

where n2n-i denotes the space of polynomials of degree 2n — 1 in one variable 
and Cw is the normalization constant so that Cw J_i w{x)dx = 1. It is known that 
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Xk > and Xk,n are zeros of the orthogonal polynomial p„ with respect to w. When 
w = Wa^p, the orthogonal polynomials are the Jacobi polynomials P^'^^ and Xk,n, 
1 < k < n, are the zeros ofP^"-'^\ We define 

Theorem 3.1. For W_i on Q, the Gaussian cubature rule of degree 2n ~ 1 is 
r n ^ ^ 

/ f{u,v)W_i{u,v)dudv = 2^^ XkXjf{uj^k,Vj^k), f^'^-jn^i, 
■'^ k=i j=i 

where means that the term for j — k is divided by 2. For Wi on il, the Gaussian 
cubature rule of degree 271 — 3 is 

n k — 1 

f{u,v)Wi{u,v)dudv = 2^^Xj^kf{uj,k,Vj,k), f e 
where Xj^k = AjAfc(x-,,„ - Xk,n)^- 



(3.3) 



(3.4) 



Proof. The product of (3.1 1 is a cubature rule on [—1, 1]^ 



(3.5) 



/n n 
f{x, y)w{x)w{y)dxdy = V V XkXjf{ 



~t~ -^j.n ) •^k.n'^ j.n) i 



which is exact for / e Ii2n-i x Ii2n~i, the space of polynomials of degree at most 
2n — 1 in either x or y variable. Applying (3.5) on the symmetric polynomials 
f{x + y, xy) and using the symmetry, we obtain 

. " / 

Cw I f{x + y,xy)w{x)w{y)dxdy = XkXjf{xk, 
■^^ k=l 3=1 

exact for polynomials / in n2n-i x ^2n~i- Under the change of variables u — x + y 
and V = xy and by (2.9|, the above cubature becomes (3.3), since n2„_i x n2„_i 
becomes ^2n-i under the mapping {x,y) i— {u,v). It is easy to see that (3.3| has 
dimn^j_j^ nodes, so that it is a Gaussian cubature rule. 

To prove (3.4), we apply the product Gaussian cubature rule (3.5) on functions 
of the form {x - yYf{x + y, xy) for / e n2„-2 x U2n-3 to get 

fc-i 

{Xj^n - Xk,nf' f{Uj,k, Vj^k)- 

k=2 3=1 

Since {x — y)'^w{x)w{y) — Wi {u, v) for u — x + y and v = xy, the above cubature 



n k—1 

I fix + y, xy){x - yfw{x)w{y)dxdy = X X '^3^k{^ 



rule becomes (3.4) under (x, y) i— > {u,v). 



□ 



By Theorem 2.1 the nodes {{xk,n + Xj^m Xk^nXj^n) ■ ^ ^ j ^ k < n} of the cuba- 



p(-5)l 



ture rule ( 3.3 1 are common zeros of the orthogonal polynomials in {Pq „ 
and the nodes {{xk^n + Xj^n, Xk,nXj^n) : 1 < j < fc < n — 1} of the cubature 



rules (3.4) are common zeros of {Pq 



n— 1 7 • ■ ■ 1 ^ n— l.n 



P. 



_]^}. Formulated in the lan- 



guage of algebraic geometry, this states, for example, that the polynomial ideal / — 



(P 



(-1) 



Pn,n^ ) has the zero-dimensional variety V = {{xk.n + Xj^n,Xk.nXj.n) 



k.n ' • ■ 



I < j <k <n}. 
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We remark that the above procedure of deriving cubature rules for Wj.i on 
[—1,1]^ can be adopted for other types of cubature rules besides Gaussian cubature 
rules. In fact, instead of starting with the product Gaussian cubature rules for 



w{x)w{y) on [—1, 1]^ as in the proof of Theorem 3.1 we can start with a product 
cubature rule of other types. For example, we can start with a quadrature rule 
of degree 2n for w that has all nodes inside [—1,1], in which case an analogue 



of Theorem 3.1 was established in [13]. We can also start with a Gauss-Lobatto 
quadrature for w to get a cubature rule that has nodes also on the two linear 
branches of the boundary of fi. 

shows that Gaussian cubature rules exist for i . An imme- 



The Theorem 



3.1 



diate question is it Gaussian cubature rules also exist for the weight function 
for 7 ^ ±i. The answer, however, is negative. 

Theorem 3.2. Forn > 1, the Gaussian cubature rules do not exist forW^i/2 -1/2 7 
z/7^±l/2. 

Proof. It was shown in [TH (10.7)] that a basis of orthogonal polynomials of degree 
n with respect to W_i/2,-i/2,7 is given explicitly by 

where < A: < n and Pn"'^' is the Jacobi polynomial of degree n. It is easy to 
see that these polynomials do not have common zeros (considering, for example, 

k = n first). Consequently, {P^. ^' {x,y) : < k < n} does not have dimn^_]^ 
common zeros for n > 1. Hence, the Gaussian cubature rules do not exist according 
to Theorem [O □ 



For what we will do in the following subsection, we make an affine change of 
variables u = 2{s — t) and v = 2s + 2i — 1, which implies that the measure becomes 
Wj{u,v)dudv = W*{s,t)dsdt, where 

(3.6) W;{s,t) ■.^2hlA'^+^w{x)w{y){{l-^lf-ty{{l + ^f-ty 

with s=\{l + x){l + y), t = i(l - x)(l - y), 

and the domain Vl becomes Vl* defined by 

VL* := {(s, t):s>0, t > 0, Vs + Vi<l}, 

which is depicted in the right figure of Figure 2. In the case of w = Wa^p the weight 
function W* becomes 

W:^^^^{s,t) = 26„,0,^4"+'^+-'+is"t^((l - ys)2 - trai + Vs)' - ty. 

Since the affine transform does not change the strength of the cubature rules, the 

Gaussian cubature rules exist for the weight functions W* 1 . Let us denote by 

* 2 

^j,k = ^j,k,n '■= i(l + + ^k,n), y*j^k — y*j,k,n •= i(l ^ •'^j,n)(l ^ Xk^n)- 

Corollary 3.3. For W^i on fl* , the Gaussian cubature rule of degree 2n — I is 

2 

« n k 

(3.7) / /(s,t)M^:,(s,0dsdt = 2^^'A,A,/(2:*,fc,y*,), f eUj^^,. 

k=l j = l 
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Figure 2. Nodes of cubature rules of degree 19 for W_ 
and W* 1 1 1 



For Wt on fl* , the Gaussian cubature rule of degree 2n — 3 is 

2 

„ n k — 1 

(3.8) / /(s,t)W^|(s,t)dsrfi = 2^^A,,fc/(4„2/*fe), /enL_3- 

The nodes of the cubature ruled of degree 19 for W_i _i _i and W* i i i 

2> 2> 2 -2'^2^^2 

are depicted in the left and the right figures of Figure 2, respectively. 

4. Minimal cubature rules 

Our goal in this section is to establish minimal cubature rules for the weight 
functions on [—1,1]^. We shall do so by several transformations of the Gaussian 
cubature rules in the previous section. 

First we recall the Sobolev theorem on invariant cubature rules. A cubature 



rule in the form of (1.1) is invariant under a finite group G if the equality is 
unchanged under / i->- af, where <jf{x) = f{xa), for all a ^ G. The Sobolev 
theorem states that if a cubature is invariant under G then it is exact for a subspace 
V of polynomials if and only if it is exact for all polynomials in V that are invariant 
under G. 



We start from cubature rules for W* in Corollary 3.3 and make a change of 
variables {s,t) The domain fl* becomes the triangle T := {{u,v) : 

u,v > 0,1 — u — V > 0} and, since the weight function W*{u'^, u^) is even in both 
u and V, we extend it by symmetry to the rhombus R, depicted in Figure 3, 

R -.^ {{u,v) : -1 < u + V < 1, ~l<u~v<l}. 

The change of variables has a Jacobian dsdt — 4\uv\dudv. We define the weight 
function on R by 

U^{u,v) -.^luvlWliu^, v^) = 2bl4'-'+'^w{x)w{y)\uv\{{l - - v^)''{{l + - v^^ , 



where u ~ i\/l + x^/l + y, v = — x^/Y 
In the case of W* a , the corresponding weight function is 



y- 



(4.1) C/a,/3,7(u,w) = 26„,^,^4 



a+0+7+ly2Q+l,, 2/3+1 



V 
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Under the change of variables (s, t) i— > (u, v) and using the symmetry, the integrals 
are related by 



(4.2) / f{s,t)dsdt^A f{u\v^)uvdudv = / f{u\v^)\uv\dudv, 

Jn* JT JR 

from which it is easy to see that satisfies Uj{u, v)dudv = 1. 

Directly from its definition, the weight function Uj is evidently centrally sym- 
metric. To state the cubature rules for Uj, we introduce the notation 9k.n by 

Xk,n ^ C0s9k,n, fc = 0, 1 , . . . , 71. 

Since w is supported on [—1, 1], the zeros of the orthogonal polynomial p„ are all 
inside [—1, 1], so that < 9k.,i < tt. 

Theorem 4.1. For U i on the rhombus R, we have the minimal cubature rule of 
degree An — 1 with diTnI{\^_^ + n nodes, 

f{u,v)U_.^{u,v)dudv^-Y^Y^\k\, 

fc=i j=i 

X ^ / cos ^ COS ± sin ^ sin , / e ^\n~-i^ 

where the innermost ^ is a summation of four terms over all possible choices of 
signs. For Ui on R, we have the minimal cubature rule of degree 4n — 3 with 
dimn|„_3 + n nodes, 



n k — 1 



„ 1 " ^ 

(4.4) / f{u,v)U.iu,v)dudv ^ 

X ^ / ^± COS ^ COS , ± sin ^ sin , / e ^\n- 
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Proof. Changing variables s = v? and t — hi (3.7) and applying (4.2), we obtain 

k 



•' ^ k=l 3 = 1 

n ^ ^ 

= ^ ^ AfcAj / {cos^ %^ cos 



2 Y>sin2%^sin2'%^ 



k=\ j=i 



for all / G n2„-i, where we have used the fact that x 



■29,, 

— sm 



j-k 



2 6i n 2 9k n 
COS COS^ . 



and 
The 



, which follows from the definition of x* and y* f, 
above cubature rule can be viewed as (4.3) applied to f{x'^,y'^). Since {/(x^,?/^) : 
/ £ ^2n-i} consists of all polynomials in 11^^-1 ^^sX are invariant under the group 
Z2 X Z2, it implies, by the Sobolev theorem, cubature rule (4.3). Since none of 
the nodes of (|3.7[) are on the edges s = or i = of 51* , the number of nodes of 



cubature rule (4.3) is exactly 



4dimn^_i 2n{n + 1) = dimHa^ 



= dimll 



2n-l 



2n 
2 ' 



which attains the lower bound in (1.4). The proof of the cubature rule (4.4) is 
similar. □ 

The nodes of the cubature rules of degree 19 for U_i _ 1 _ 1 are depicted in the 

2 ■ 2 ■ 2 

right figure of Figure 3. 

As a final change of variables, we rotate the rhombus by 45° to the square 
[— 1, 1]2. This amounts to a change of variables u = {x + y)/2 and v = {x — y)/2. 
The measure under this change of variables become J7^(ti, v)dudv = yV^(a;, y)dxdy, 

Wj{x, y) =6X4'^w(cos(6i - 0))w(cos(6i + <j)))\x^ ~y^\[l~ x^V{l - y^)\ 
where a; = cos 6*, y = cos(f), (x, y) e [— 1, l]'^. 

A simple computation shows that this is precisely the weight function defined in 
(2.12). In the case of Ua^p^-y, the corresponding weight becomes Wa^p,-/ defined by 

n'^,p,j{x,y) = 6aA74^|x + 2/P"+i|x - y\'^+\l - x^^l - y^)\ 



which is exactly (2.13). Since the strength of the cubature rules do not change 
under the affine change of variables, we then have minimal cubature formulas for 
VVa,;3,7- To state this cubature explicitly, let us define 



(4.5) 



Sj,fc 



cos ■ 



and tj^k '■— cos 



9j.n+0k^' 



where Ok.yi is again the angular argument of the zeros Xk.n = cosf?fc.„ of p„. 

Theorem 4.2. For yV_ i on [—1, 1]^, we have the minimal cubature rule of degree 
An ~ 1 with dimn2„_]^ + n nodes, 

fix, y)W^i{x, y)dxdy ^ A^A, [/(sj,fe, t,, fc) + /(tj,fe, Sj,fc) 
.-1.1]^ 2^1 7^1 

+ f{-Sj,k, -tj,k) + f{-tj,k, -Sj.k)] ■ 
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For Wi on [—1,1]^, we have the minimal cubature rule of degree 4n — 3 with 
diinn2,j_3 + n nodes, 

n k—l 

f{x,y)Wi{x,y)dxdy =^^^1^^.^ [f{sj,k,tj.k) + f{tj,k,s,,k) 

+ fi-S],k, -tj,k) + f{-tj,k,-Sj,k)] , 

\2 



where Aj.fc ~ XjXk{cost 



cos Uk ; 



1 _ 1 _i_ 1 , these cubature 

2 ' 2 2 



In the case of the product Jacobi weight function W, 
rules were constructed in and, more recently, in [5] via a completely different 
method. In all other cases these cubature rules are new. The nodes of the cubature 
rule of degree 35 for the weight function 



are depicted in the left and right figures in Figure 4, respectively. The influence 





Figure 4. Nodes of minimal cubature rules of degree 35 for 



and W, 



0,0,- 



of the part |a;^ — y^j in the weight function Wq q _i is clearly visible in comparing 
with the cubature rules for W_ i _ i _ i . 

2' 2' 2 

By the relation (2.12 1 and the integral relation (2.14), we could arrive at cubature 
rules (4.6 ) and (4.71 from those in (3.3 ) and (3.4 1 by the mapping (x, y) i— (2a;y, x^-l- 
2/2 — 1)^ bypassing some of the middle steps. Our presentation, on the other hand, 
is more intuitive and provides, hopefully, a better explanation of the connection 
between the Gaussian cubature rules for and the minimal cubature rules for 



We can also give a proof of Theorem 4.2 based on Theorem |2.3| by consider- 
ing the common zeros of corresponding orthogonal polynomials, although a direct 
computation of the cubature weights will not be easy. Recalling the orthogonal 



polynomials iQ^^Xn defined in (2.151, the following corollary is an immediate con- 
sequence of Theorem |2.3| 
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Corollary 4.3. The nodes of the minimal cuhature rule (4.6) a re th e common zeros 
of orthogonal polynomials {iQ'^'2n '■ < k < n} in Proposition 



2.7 And the nodes 



of the minimal cubature rule (4.7) are the common zeros of orthogonal polynomials 

{iQt:2;L--o<k<n-i}. 



The relation (2.15) shows that the nodes of the minimal cubature rule (4.6) and 
the nodes of the Gaussian cubature rule (3.3) are related by a simple formula: if 
{s,t) is a node of the former, then {2st, s^ + t^ — 1) is a node of the latter; fur- 
thermore, the nodes {sj^k,tj,k), {tj,k, Sj^k), {-S:i,k, -tj^k), {-tj,k, -Sj,k) of the former 
correspond to the same node {2sj^ktj,k, ■^j fc + fc ~ 1) latter. This can also 

be verified directly by elementary trigonometric identities. 

It should be pointed out that the Theorem |3 . 2| shows that the above construction 
does not work for the weight functions 



(x,y) = {l-x'ril-y'r 



when 7 7^ ±1/2. We cannot, however, conclude that the cubature rules of degree 
An—1 that attain the lower bound (1.4 1 do not exist for these product Gegenbauer 



weight functions. In fact, examining the proof carefully shows that the procedures 
that we adopted could be reversed only if the cubature rules for yV_ \ satisfy 
certain properties. What we can conclude is then the following: If a cubature rule 
of degree 4n — 1 that attains the lower bound (1.4) exists for W_i -i 7: then either 
it is not invariant under the symmetry with respect to the diagonals y — x and 
y — —X of the rectangle [—1, 1]^ or some of its nodes are on these diagonals. 

Finally, our procedure of deriving cubature rules for Wj-i on [—1,1]^ can be 
adopted for other type of cubature rules, such as cubature rules of even degree or 
Gauss-Lobatto type cubature rules, see the remark at the end of Subsection 3.1. 
In particular, if we start with a Gaussian-Lobatto quadrature for w, which has 
additional nodes at —1 and 1, then the resulted cubature rule for W-y will have 
nodes on the diagonals of [—1,1]^. Since they do not seem to have other features, 
we shall not pursue them further. 



5. Lagrange interpolation and Gaussian cubature rules 

Cubature rules are closely related to Lagrange interpolation polynomials, as 
stated in Section 2. In this section we consider Lagrange interpolation polynomials 
based on the zeros of the Gaussian cubature rules constructed in the Section 3. 

The Lagrange interpolation polynomial based on the Gaussian cubature rule in 
Theorem |2.1| is given in Theorem |2.2| A more direct construction can be given 
however as follows. 

Let {a;fe,„ : 1 < A: < n} be the zeros of the orthogonal polynomial p„ of degree 



n with respect to w on [—1,1], as in (3.1 ). The Lagrange interpolation polynomial 

Vn{x) 

Pni^k 



(of one variable) of degree < n based on these points is 

n 

(5.1) lnf{x)^^y2^f{xk.n)lk{x), lk(x) := —r- 

^1 P'niXk,n){x - Xk,n) 
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Theorem 5.1. The unique Lagrange interpolation polynomial of degree n—1 based 
on the nodes of the Gaussian cubature rule (3.3) is given by 

n k ^ 

(5.2) L„/(m,u) = ^^ f{uj^k,VjM)lj,k{u,v), 

k=l J=l 

withlj^k{u,v) ■.^lj{x)lk{y) +lj{y)lk{x), u = x + y, v = xy. 

And the unique Lagrange interpolation polynomial of degree n—2 based on the nodes 
of the Gaussian cubature rule (3.4) is given by 

n k — 1 

(5.3) Lnf{u, V) = ^Y1 fi'^hk,V;j,k)lj.k{u, v), 

■ , / X h(x)lk{y) - l]{y)lk{x) 

with L k(u,v) := , u = x + y,v — xy. 

x-y 

Proof. A quick computation shows that if < j < fc < and < p < q < n, then 

^jyki^p^qj ^p^q) ^ j i,*^ p ^n)^ ki,*^ q .n) ~^ ^ j {.-^ q ^ri)^ ki,-^ p ^j^p^k^q ~t~ ^k^p^j,q ^j^p^k^q- 

liO<j — k<n and < p < q < n, then 

,k{'^p.q^ ^p.q) — '^lk{-^p,n)lk{'^q,n^ — ^^j^p^k.qi 



which proves (5.2). The proof of (5.3) is similar. 



□ 



The exphcit formulas of lj,k can also be obtained from Theorem 12.21 In fact, as 

iTTm (±-) 

shown in [551 Theorem 3.1], the reproducing kernel Kn ^ (•, •) = Kn{W^^] •, •) can 
be expressed in terms of the reproducing kernel 



kn{x,y) = kn{w;x,y) := ^Pk{x)pk{y) 



k=0 

of one variable, where pk are orthonormal polynomials with respect to w. Set 

u := {ui,U2) ^ {xi + X2,xiX2) and u (ui, U2) = (yi + 2/2, 2/12/2)- 
The reproducing kernel Kn ^ (•, •) ^or W_i is given by 

(5.4) Kn ^{u,v) = -[kn{xi,yi)kn{x2,y2) + kn{x2,yi)kn{xi,y2)], 

(-) 

and the reproducing kernel Kn {■, •) for Wi is given by 

fK K\ T^^^h \ - kn+l{xi,yi)kn+l{x2,y2) " fc„+l (2:2 , 2/1 ) fcn+1 (a^l , 2/2 ) 

[0-0) J\n [U,V — — -- r . 

2(xi - X2){yi - 2/2) 

As an application of the explicit expression, we can estimate the uniform norm of 
the interpolation operator, often called the Lebesgue constant. For the interpolation 
polynomial /„/ in (5.1), the Lebesgue constant |l/n||c[-i,i] satisfies 

n 

ll-^n||c[-i,il = max V|^fc(a;)|. 
Corollary 5.2. The Lebesgue constant for Lnf in (15. 2| satisfies 



(5.6) 



Ijinlloo < (||^n||c[-l,l])' 
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Proof. A standard argument shows that the Lebesgue constant for £„/ is given by 

n ^ 



^mioo-.max \lj^k{u,v) 



(M,ti)60 



/c=l 3 = 1 



Since tj^k{u,v) — £k,jiu,v) by (5.2), a moment of reflection shows that 



J2 E ^) H E E 1^^- (^)'^ (j') + (2/) I < 2 E 1^^- (^) I E 1^^- (y) I ' 

fe = l J = l fc=li = l j = l 3 = 1 



from which the estimate (5.61 follows immediately. 



□ 



Denote by L^"^/ the La grange interpolation polynomial based on the nodes of 
the Gaussian cubature rule of degree 2rt — 1 for p - 

Corollary 5.3. Let a,l3 > —1. The Lebesgue constant of L"'^ f satisfies 

_ max{a,/3}>-l/2, 
I log n, maxja, p| < —1/2. 

Proof. This follows from the previous corollary and the classical result on the La- 
grange interpolation polynomials at the zeros of Jacobi polynomials [16] . □ 

6. Lagrange interpolation and minimal cubature rules 

The relation between a minimal cubature rule and the Lagrange interpolation 
polynomial based on its nodes is stated in Theorem |2.5[ In this section we discuss 
the Lagrange interpolation polynomials based on the nodes of the minimal cubature 
rules of degree An — 1 in Section 4. In order to derive explicit formulas and discuss 
the Lebesgue constants, we shall limit our discussion to ^ - j ' which we renamed 
as Wa,/3 at (2.22). An analogue discussion can be carried out for ^ i . 



6.1. Construction of the interpolation polynomial. Let X„ denote the set 
of nodes of the cubature formula Wa.fj- The Lagrange interpolation polynomial 
based on X„ is given in Theorem : 



polynomial P, 



(",/9) 



2.5 



in which „ — x^"'^^ are the zeros of Jacobi 



The subspace IIj^ in (2.2) now takes the form 



TlL-i U span{2Q^"2?/ : < fc < n - 1}. 

The interpolation polynomial in IIj^ is given in Theorem |2.5| in terms of a kernel 
of Iljn defined by 



(6.1) 



fe=0 



where bk,n are certain positive numbers, JC2n-i and 2Qk.2n are given explicitly in 



(2.23) and Proposition (2.7) 



Although the cubature rule (4.6) of degree 4n — 1 for W_i can be deduced 



from the Gaussian cubature rule (3.3) for W_i, this deduction does not extend to 



interpolation polynomials, since each node of the cubature rule (3.3) corresponds 



to four nodes of the cubature rule (4.6 ). We have to work with the explicit formula 



given in Theorem |2.5[ which we determine explicitly in the following theorem. 
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Theorem 6.1. Let Xk,n — cos9k,n, ^ l£ k < n, denote the zeros of the Jacobi 



Set 



polynomial P^^'^^ and let Sj_k '■— cos ^j-" ^^''■^ and t 
Then the Lagrange interpolation C"'^ f in H* is given by 

n k 

(6.2) c:/f{x,y) - EE [/ i^i) 4!i(^'y) + / (^5-5) 45(^-2^) 

fc=i j=i 

+/(x55)4^(x,z/)+/(x2)4:^(->^)^ 

where the fundamental interpolation polynomials fj'], are given by 



(6.3) 



2n 



((a;,y),xg) , 



in which | in the right hand side needs to be replaced by j when j — k, and 
(6.4) ^L(x,y)=^,"f_,(.,t)+,l + " + '^ + "^(^'^)^^^ 



l + a + 



J^d^::;\xi-xi){yi-yi) 



K 



Q + 1,/9+1 
n-1 



{s,t)-K 



Q + 1,^+1 



(s,t) 



n{l + a + fi + n) 
"(l + a + /3 + 2n)2 



n-2 

2(5n-l,2n(a;)2(3 n— 1 ,2n 



w/iere s = {2x iX2,xl + - 1), t = {2yiy2, yt + ? /| - 1), IC2„_i{-, ■) an-rf d^,}? 
given in (2.23) and m _R'"''^(-,-) is given in (2.24). 



are 



Proof. The formulas (6.2 1 and (6.31 are exactly those given in Theorem 2.5 spe- 
cialized to the Jacobi case. It remains to establish the formula of (6.4), for which 
we need to determine the constants bk n in (|6.1| . 



Throughout this proof, we write Qk.2nix,y) = 2Q'k'f\^Ty)- By the explicit 



formula of Qk,2n in Proposition |2.7[ it is easy to verify that 
(6.5) Qman i^fl) = 1»,I3\J'^- ^j\A 



and furthermore, since Qm,2n is symmetric in its variables, 



(6.6) 



Wm.2n \ ^j^k } ' 



Let us denote by C„[/] the minimal cubature rule, that is, 

k 



CM] :='EE'^ 

k=i j=i 



- -k --J 



(",/3);^(" 



By (2.4) and the fact that are the fundamental interpolation polynomials, we 
obtain 



'^2n y^j,kT-^j',k' 

which implies immediately that 



= 2 



(6.7) 



c 



2« (Xj^i, •) Ql,2n - Qu2n [^f]^ 
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On the other hand, using the formula of /C2n('j ') (6-1) shows that 



n— i 

^2n i^^f.h ') *3'.2n = C„ -f^2n-l ("'^j^fc' ') Qi.Znj + bm,nCn [Qm,2nQl,2n 



m=0 



Since the cubature rule is of degree An — 1 and Qi.2n is an orthogonal polynomial 
of degree 2n, 



Cn 



K2rf-i (xg, •) QL2n\ = 2^^^ J ^ ^ K^f_, [^^^ y) QiMy)Wc.Ay)dy 



= 0. 



Furthermore, since Qm,2n is symmetric in its variables, it follows from (6.6) that 

n n 

Cn [Qm,2nQl,2n] — -^fc ^'"^i QirL,2n (-""Ifc) Ql.2n j 



Recall the definition of h,n defined in (2.20). By (6.6), the explicit formulas of Qk,2n 
and the Gaussian quadrature (2.19), it follows that 

Cn [Qm,2nQl,2n] = ^^l^phn-lKuSlm, < l,m < n - 1. 

Putting these formulas together, we have shown that 



r 



^2n i^i L'l QL2n ^ "^la.B^n-l 



Comparing with ( |6.7| , it follows readily that &j ^ = 27^ phn-ihm- Recalling that 



i,/3+i/(\/2ca,^), applying Lemma 



2.6 



gives 



^0,n — ■ • • — bn-2,n — Z — ; 

1 + a 



1 + a + /3 + n 
+ a + (3 + 2n' 



and 6„_i,„ = 60, n • 



The final step in verifying (6.4) uses the fact that 

n-l 



k=0 



which can be verified using the explicit formulas of the quantities involved and the 
elementary trigonometric identity 



, x'^ + y"^ — 1 — cos( 



(6.8) 2xy = cos(6' -(/)) + cos{0 - 
see also Section 4 of [22]. This completes the proof 



I cos(6' + ( 



□ 



The above theorem gives a compact formula for the Lagrange interpolation poly- 
nomial based on the nodes of the minimal cubature rule with respect to Wa,/3- In 
the case of a — (3 = —1/2, the interpolation polynomials were introduced in [20] 
and they were studied numerically in The explicit formulas given in [20], how- 
ever, takes a different form since the set of nodes were not divided into the four 
subsets as in (6.2) and a completely different formula for /C2n('i ') used. 
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6.2. Lebesgue constants of the interpolation operator. The Lebesgue con- 
stant of the interpolation operator C^'^ is its operator norm ||oo- Since 

||/:^:'^/||oo < lir^'^^ooll/lloo, V/eC[-l,l]2, 
the Lebesgue constant determines the convergence behavior of 



Lemma 6.2. The Lebesgue constant of Cf^'^ f in Proposition 6.1 satisfies 

k 



(6.9) ii^n^iu = LerMpEE'^.- 



k=l J=l 



/Co 



2n \ -^^^j.k 



(2) 

2ri \.^J^,,fc 



max 

a;e[-l,l] 



(--42) 



fe=l 3 = 1 



Proof. RecaUing (4.5) and the definition of ^, it follows easily from the symmetry 

that ||/3"'^||oo is bounded above by 4 times of the quantity in the last expression 
and it is at least as big as the same quantity. □ 

In order to deduce the order of the Lebesgue constant, we need to estimate, by 
the explicit formula at (6.4), several sums. We first deal with the easiest sum to be 
estimated. Let c denote a generic constant whose value may vary from line to line. 

Lemma 6.3. For a, (5 > —1 and x ^ [—1,1], 



Ac 



k=ij=i 



2Vn-l,2n 



(a;)2g„-i,2„ (a;,xg)| <cn2--{«^«. 



Proof. We will need several well known estimates for the Jacobi polynomials and 
related quantities, all can be found in [16 . First we need 

(6.10) b^«(x)|<c(^^^+n-l)-^"+^/^'/^(V^+^+n-l)-('^+^/^^/^ 

for a; € [—1,1]. Using the fact that cos^ 9 — cos^ 4> — sm{9 — (j>) s'm{d ^ 
from the explicit expression of 2Qn-i.nix) that 

|2Q„-i,2„(x)| <c max iVl -a;2pi":«(x) I <C7^2 xe[-l,l] 

—l<x<l 

Furthermore, we need the estimates 



it follows 



.11) 



.12) 



A 



k.n 



~^Wa,l3{xk,n)\ 1 - 



Pn-l'{^k,n) ^ [Wa,0{xk,n)] ^ (1 - Xk,n) 



From (6.12), it is not difficult to see (using (2.21), for example) that 
Consequently, since 2Qn-i,2n is symmetric in its variables, we see that 



as the sum is easily seen to be bounded upon using (6.11). 



□ 
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The other sums of ||>C^'^||oo cannot be deduced form the Lebesgue constant for 
the interpolatfon polynomial of one variable, as we did in Corollary |5.2| since there 



are four remaining sums by (6.4 1, and only one of them, the first one, is related 
directly to the fundamental interpolation polynomials of one variable. We can, 
however, reduce the proof to the estimate of several kernels in one variable. Let us 
define, for i, j > 0, 



:= (1 - + - y)4(l 

Lemma 6.4. Let a, 13 > —1/2 andi.j > 0. Then 
(6.13) |fci"''^)''''J'(cos6',cos0)| 



< c 



(sinfsinf 01 + n-'^)-°'~i {cos ^ cos ^ 

|0-(/)|+n-i 



While (6.11 ) and (6.11 1 are classical, (^6.13^ is stated recently in 



and its proof follows from an estimate in The restriction a,f3> 
a,/3 > — 1 comes from the method used in 4 . For i,i > 0, let 



Lemma 5.3] 
-| instead of 



k=l 

We will need the following result for our estimate of ||£"''^||oo. 
Lemma 6.5. Let a, l3 > -1/2. Fori,j>0, 

-^max{a,/3}+i^ max{a, /?} > -1/2, 



(6.14) 



max A(;'^)(x) = 0(l) 




ma,x{a, /3} = —1/2. 



Proof. We can assume x G [0, 1] and write x — cos 6. We consider a > —1/2, the 
case a = — 1/2 is easier. Fix m such that Xm,n is (one of) the closest zero to x. 
Then l<rn<n/2 + l. We only consider the sum in A^''^^ for 1 < A: < 2n/3, the 
remaining part is easier since for 2n/3 < k < n, |0— 0fc,„| ^1. If fc = rn — 1, m, m+1, 
then by and ([oTl), 



A 



< 



(singfc,„)°+2 
(sin^ ^ +n-2)"+3 



< cn 



Using the fact that |0 — i 



E 



|fe-m|>l 
l<fc<2n/3 



9,n - 0|, we have by ( [6l3| and ( [6TT| ) 



|fc-r?j.|>l 
l<fe<2ri/3 



to|(/cto 



The last sum can be shown to be bounded by dividing it into three sums over 
1 < A; < to/2, to/2 < k < 2m and m < k < 2n/3, respectively. Such estimates are 
rather standard affairs, we leave the details to the interested readers. □ 



For i — j — 0, the estimate (6.141 gives the order of the Lebesgue constant for the 
interpolation polynomials based on the zeros of Jacobi polynomials in one variable. 
The classical proof in [Hj, however, does not apply to the case of ^ (0,0), 
since A^"„''^fci"'^'''*'"' (a;, Xfe^„) does not always vanish at xi,n when I ^ k. 

We are now ready to prove our result on the Lebesgue constant of C'^'^f- 
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YUAN XU 



Theorem 6.6. Let a,/3 > —1/2. The Lehesgue constant of the Lagrange interpo- 
lation polynomial C"'^ f based on the nodes of the minimal cubature rule of degree 
An — I for yVa.p satisfies 



(6.15) 



i^^:^''iioo-o(i) 



^2 inax{a,/3}4-l 
(l0g?l)2, 



max{Q!, (3} > -1/2, 
maxja, /?} = —1/2. 



Proof. Let Xk^n = cos 0k,n be the zeros of the Jacobi polynomial p-h . We estimate 
ll'C^'^lloo in (6.9) by setting xi = cos and X2 = cos and taking the 

maximum over < 0i, 6*2 < tt. It follows that 2a;ia:2 — cos 0i + cos 02 and x\-\-x\ — 



1 = cos 01 cos 02, and furthermore, 
Xi — X2 = ^(1 — cos0i)(l — cos 02) 



and 



Xi 



+ X2^ + COS0i)(l + COS 02). 



Hence, recalling (4.5), it follows from (6.4), (2.24) and Lemma 6.3 that 



O 1) max W XiXk 



J 



0,1/ 



J, 



1,1/ 



1,02) 1 +0(1)712 



where J^'-' are defined by 

J]i{e^,02) (cos 01, cos 0,)fci"'^)'''^"(cos 02, cos 6k) 

+ fci"'^)'*'^' (cos 01 , COS 0fc)fci"''3) '^J' (cos 02 , COS 0j ) . 

Hence, as in the proof of Corollary |5.2[ we can reduce the estimate to the product 



of K^n '\ SO that the desired result follows from (^6.14h. 



□ 



In the case oi a — (3 — —1/2, the order of the Lebesgue constant was determined 
in [3] based on the explicit expression in [5D]. In all other cases, the estimate (6.15) 
is new. One interesting question is if the result can be extended to the case of 
max{a,/3} < ~\- We expect that it can be and, furthermore, we believe that the 
order is ||/:""^||oo = 0(l)(logn)2 for max{a,/3} < -\. 
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